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ABSTRACT 

The formation of dust with temperature-dependent non-grey opacity is considered in a se- 
ries of self-consistent model atmospheres at different phases of an O-rich Mira variable of 
mass 1.2 Mq. Photometric and interferometric properties of these models are predicted under 
different physical assumptions regarding the dust formation. The iron content of the initial 
silicate that forms and the availability of grain nuclei are found to be critical parameters that 
affect the observable properties. In particular, parameters were found where dust would form 
at 2-3 times the average continuum photospheric radius. This work provides a consistent phys- 
ical explanation for the larger apparent size of Mira variables at wavelengths shorter than 1 /xm 
than that predicted by fundamental-mode pulsation models. 

Key words: techniques: interferometric - stars: variables: Miras - stars: AGB and post-AGB 



1 INTRODUCTION 

Dust shells around Mira variables often dominate the mid- and 
far-infrared regions of their spectra and strongly influence the ob- 
served brightness distributions of their stellar and circumstellar ma- 
terial. In contrast to C-type Miras for which substantial theoretical 
progress has been made in recent years (e.g. Helling et al. 200i 



Schirrmacher et al. 2003; Hofner et al. 2003; Gautschv-Loidl et al. 
2004), for oxygen-rich M-type Miras the physical conditions and 
the geometric region of formation of dust grains is poorly under- 
stood. 

This situation is first of all due to the much more complex 
physics of formation of typical dust grains, which may contain 
metal oxides such as corundum (AI2O3), silicates or even solid iron 
lGaill2 003). These heterogenous grains form around nuclei that are 
in general c omposed of different compounds again, such as Ti02 
Jjeong et aj2003h . Furthermore, in order to understand the occur- 
rence of dust in an M-type Mira, the formation of grains has to 
be studied within the physical environment of the time-dependent 
dynamic atmosphere of the pulsating star. Interferometric observa- 
tions indicate that the inner radii of dust shells may be as small 
as 2 to 4 stellar radii where the stellar radius refers to the po- 
sition of the continuum-forming layers (e 
Danchi & Be"sterl Il995t Lopez et all \\w\ 
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Irela nd et alj2004dL l2005V This is is well within the upper regions 



the ste llar spectr um are form ed (cf. I Hofmann et all 1 9981: iTeiet alJ 
12003 JtilOnnakj2004IPerrin et al 120041: lOhnaka et all2005h 

So far existing dynamic model at mospheres of M-type Mira 
stars either are dust-free (e.g. IHofmann et allll998t IWoitke et alJ 
Il999t HofneretalJ 2003), or they use grey, composition in depen- 
dent dust opacities (e.g. IWinters et alj|200ot Ijeong etalll2003l) . 
iBedding et alJ 1200 ll) studie d schematically the co nditions of dust 
formations in the models of IHofmann et alJ 1 19981) and concluded 
that corundum and silicate grains are likely to occur, depending 
on the phase of pulsation, in the upper layers of the star's atmo- 
sphere. However, only a small number of observable features of a 
small number of models were considered, with a very simple opac- 
ity treatment. 

In this paper we investigate the formation of dust in a series of 
dynamic models of a Mira variable, and discuss which observable 
effects would result from different physical assumptions regarding 
the dust formation. In particular, we focus on effects of non-grey 
opacities of silicate dust with a variable Fe content and the effect 
of grain size on observable properties. As hetrogenous O-rich dust 
formation is complex, our approximations are designed to be most 
valid for only the innermost dust that has a significant influence on 
observable properties. We aim to parametrize uncertainties in this 
first significant dust formation with the smallest number of physical 
free parameters. 



of the stellar atmosphere in which typical molecular features of 
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2 SIMPLE DUST MODELS 

In this section, we will describe the procedure we used to calculate 
absorption and scattering coefficients for silicate dust. We will then 
discuss the degree of validity of this procedure for the physics ap- 
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plicable to Mira atmospheres, the way in which corundum dust is 
added to the models and finally the definition of the dust types used 
in the models. 

The properties of a given dust type are considered to be wholly 
determined by the dust temperature Td and the gas pressure P. Td 
is found by solving the equation of radiative equilibrium: 



K, u J u dv 



K v B v dv 



(1) 



where J„ is the mean intensity, B v = B v (Td) is the Planck 
function, and k v = n v (Td, P) is the dust absorption coefficient in 
the layer with local dust temperature Td and local gas pressure P. 
The mean intensity is found iteratively in the modelling procedure 
described in Section [3] The gas temperature T g does not enter in 
to this equation because the dust energy exchange with the gas is 
negligible compared to the dust energ y exchange with the radiation 
field (e.g. see Appendix D of Gauaer et al. 1990). 

2.1 Silicate Dust Properties 

We assume that the silicate that condenses first is olivine, 
Mg2a;Fe2-2a:Si04, and use a s implified vers ion of the non- 
equilibrium treatment of lGail & Sedl mav3il999l) to determine the 
condensation fraction of Si, /si, and the composition parameter x. 
fsi is assumed to follow the following simple relationship: 



fsi 





0.25 

1240+46 log P-T d 
100 



T d > 1240 + 46 log P 
T d < 1215 + 46 log P 
Otherwise 



(2) 



The stability limit of silicate in this equation c omes from a lin- 
ear fit to the stability limit for forsterite shown in lGaill l2003h for 
solar metallicity. The upper limit for Si condensation of 0.25 was 
chosen somewhat arbitrarily, so that there was not unrealistic grain 
growth in outer atmospheric layers, where low gas pressures mean 
that the silicate condensation fraction is far from equilibrium. The 
dependence of the olivine stability temperature on x is neglected in 
Equation[2] However, its inclusion would not cause any significant 
changes to this paper, as the lowest value of x considered is 0.775, 
where the difference in stability temperature of the olivine from 
pure forsterite is only 30 K. More sophistication was not warranted 
in this relationship because the true equilibrium condensation frac- 
tion is a function of both Td and T g , and the system is not in general 
close to chemical equilibrium. 

In practice, the Si condensation fraction fsi is determined by 
the functions x(fsi) and Td(x), given that the opacity of silicate 
around 1 is a very strong function of the composition param- 
eter x. In equilibrium condensation calculations, x « 1 for the 
initial condensate as it is energetically more favorable in the bulk 
solid to have Mg rather then Fe anions. However, the pulsating at- 
mosphere of a Mira variable is certainly not one where chemical 
equilibrium considerations apply. To find x(fsi), we consider x to 
be defined by the conditions during rapid grain growth , approxi- 
mate K p = in Equations (55-57) of iGail & Sedlmavrl 1 19991) (a 
reasonable approximation as seen by the near-constant x(Td) func- 
tions in Figure 3 of that paper) and assume that the grain has the 
same composition x throughout. This gives: 
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(4) 



Here pm s and pf c are the partial pressures of Mg and Fe left 
in the gas mixture, and msio and m.M g are the masses of an SiO 
molecule and an Mg atom respectively. Equation |4] is solved for 
x, which is limited to b e less than or equa l to 1. The elemental 
abundances e come from lPalme etail ll98ll) . for consistency with 
the modelling code, a is the ratio of the exchange coefficient for 
Mg and Fe ions to the sticking coefficient for SiO molecules. 

For a greater than or equal to 1, the assumption that the grain 
has the same composition x throughout is justified if the diffusion 
timescale for Fe and Mg anions is small compared to the timescale 
for grain growth. From the discussion in Section 5.2 o r lGaiill2003» . 
this is the case for olivine at dust temperatures above about 900 K 
and grain sizes less than about 0.1/^m. However, for small val- 
ues of q, the value of x calculated using Equation |4| will only be 
the composition of newly deposited material on a growing grain, 
and not the composition of the whole grain. During grain destruc- 
tion, Equation|4|does not accurately describe the grain composition, 
but the general property predicted by this equation that x will in- 
crease as the grain evaporates is consistent with experimental stud- 
ies lOzawa & Nagaharal200d) . 

In practice, Equation|4|is a good approximation for x during 
initial condensation or a phase of rapid grain growth if a is greater 
than about 1.0. In other cases, it represents a simple parameter- 
ization of x that includes an appropriate decrease in x as Mg is 

depleted from the gas phase. 

In forming Equation[4]from Equation 57 of IGail & Sedlmavrl 

ll999h . it was assumed that Q MgHMg = QFefFe, the sam e assump- 
tion used to create Figure 3 of lGail & Sedlmavn (j/)99 ). In this re- 
lationship, the «M g and «f c parameters are sticking coefficients 
and the v parameters are thermal velocities. This is a necessary as- 
sumption given the lack of experimental data for olivine sticking 
coefficients with intermediate values of x. 

Absorption and scattering coefficients are scaled as in the 
Rayleigh limit, where at a constant condensation fraction the ab- 
sorption coefficient is independent of grain radius, and the scatter- 
ing coefficient is proportional to the cube of the grain radius, a. The 
grain radius relates to the Si condensation fraction fsi and the num- 
ber of available nuclei per hydrogen atom iV nuc by the formula: 



3e Si Vo 

47T/ S iA^ nL 



\l/3 



(5) 



Here es; is the abundance of Si at solar metallicity, and Vb 
is the volume of the monomer, taken to be 7.5 x 10 -23 cm 3 . 
The optical constants of the olivine with x = 0.5 are taken from 
iDorschne r et alJ fl995j ). and those for forsterite with x = 1.0 taken 
from ljager et alJ 120031) . These constants are used to derive K0.5, 
fjo.5, Ki.o and cti.o in the Rayleigh limit for 30 nm grains, with the 
subscripts 0.5 and 1.0 representing x for the olivine. The absorption 
and scattering coefficients for the olivine are then approximated by 
a simple interpolation: 



Col 



/ Si ((2-2x)Ko.5 + (2x-l)Kl.o) 



(6) 



2.6 x 10" 

iVnuc 



-/ s 2 i((2 - 2x)a . 5 + (2x - l)<n. ) (7) 



2.2 Dust Formation in a Typical Mira Atmosphere 

In the above parameterisation of absorption and scattering coeffi- 
cients, four important aspects of grain growth and destruction were 
not discussed: dust nucleation, the approximation of using a sin- 
gle grain radius a, the timescale of grain growth and destruction 
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compared to the pulsation timescale and the possibilities for fur- 
ther grain growth beyond /si = 0.25. We will discuss these one at 
a time, with specific reference to the physics of a Mira atmosphere. 

The first stage of dust formation is the nucleation of seeds in a 
super-saturated gas. The number of seed nuclei available for grain 
growth determi nes the grain size. In the detailed nucleation dis- 
cussion of Jeon g et aT] 1200 3l) . the seed nuclei are found to form 
from Ti02- The nucleation rate is a complicated function of tem- 
perature and pressure: at a typical Mira atmosphere dust formation 
pressure of 0.1 dyn/cm~ 2 , the Ti02 nucleation rate increases from 
10" 29 to 10" 17 per H atom per second as temperature decreases 
from 1200 to 1150K. This relationship has significa nt uncertain- 
ties, as it is based on ab initio calculations from ljeong et alJ J2000h 
which predicts an erroneous zero-point energy of Ti02 by 0.73 eV 
or 5%, much larger than experimental errors (compare Table 2 of 
that paper with Baldu ccTet alJ ll985l) ). If one assumes all energies 
are erroneous by this relative amount, then errors in nucleation tem- 
peratures of the order of 50 K would be expected. Given these diffi- 
culties, we will consider that grain nucleation is complete when the 
growth of dominant dust components begins, and leave the number 
of grain nuclei per H atom N nuc a free parameter in our models. 

By separating nucleation and growth in this artificial way, we 
make the approximation that all dust grains in a single atmospheric 
layer are the same size. The effect of this dust on absorption and 
scattering is only equivalent to a distribution of particle sizes in the 
Rayleigh limit when the grain radius a « A, a condition always 
met for the dust types and wavelengths considered here. This ap- 
proximation is also more desirable than fixing a single particle size 
distribution at all atmospheric layers, as that approach doesn't in- 
clude the increase in scattering coefficients with increasing grain 
radius. 

As material cooling behind a shock in a Mira atmosphere 
moves outwards, the dust condensation fraction given by Equa- 
tion [2] rapidly increases as the central luminosity decreases near 
minimum and the radiation field becomes increasingly geometri- 
cally diluted. For the dust treatment here (assuming fsi is a func- 
tion of the instantaneous variables Td and P), the grain growth 
m ust be rapid compared t o the dynamical timescale. According 
to lGail & Sedlmavn 1 19991) . the limiting factor in the growth rate 
of olivine is the deposition of an SiO molecule on the grain sur- 
face. Considering a super-saturated gas mixture at HOOK with 
an assumed sticking coefficient for SiO of 0.1, the rate of grain 
growth is greater than lOnm in 0.1 cycles (3 x 10 6 s) whenever 
log(P) > —2.22. Therefore, we shall be cautious in interpreting 
any dynamical grain growth and destruction effects that occur at 
log(P) < -2.22. 

At dust temperatures where olivine has mostly condensed, 
there are a variety of possibilities for further dust condensation, 
including quartz, and solid iron. In a non-equilibrium scenario, fur- 
ther model condensation will depend on the initial choice of a, as 
well as other unknown parameters. Note that at the near-maximum 
model phases, fsi = 0.25 is never reached in the outer layers (near 
the arbitrary 5 R p model surface: see Section[3j. and the low gas 
pressures would mean that any growth would be slow compared to 
the dynamical timescale. 

2.3 Addition of Corundum 

It ha s been suggested (e.g. Maldon feTal] l2005t lEgan & Sloanl 
l200ll) that corundum (AI2O3) plays a significant role in radiative 
transfer processes in O-rich AGB stars. This is due to the clear 
need for opacity in the 12-15 /im range when fitting models to ob- 



served spectral energy distributions. However, physical models of 
Mira atmospheres with realistic gas densities (e.g. iBedding et alJ 
1200 ll) have previously shown that the opacity of corundum and the 
abundance of Al are not high enough to have a significant effect 
on observed spectra. Therefore, we will re-examine the effects of 
corundum in this paper. 

The difficulty in modelling the growth of corundum grains is 
that nucleation and growth will certainly be simultaneous, due to 
the high corundum stability temperature. In this paper we are not 
concerned about the details of this process, but attempt to allow 
corundum to have the maximum plausible effect on the radiative 
transfer. 

We have taken the stability limit for corundum at solar metal- 
licity from Gail 1 2003), fitting a function linear in log P for —4 < 
log P < 0. For corundum dust only, the Al condensation fraction 
f A i is given by: 

!0 T d > 1600 + 60 log P 

1.0 T d < 1500 + 60 log P (8) 

1600+60 log P-T d Qtherwise 

The dust temperature Td is found using Equation 1, wit h the 
optical constants of corundum taken from Koik e et all 1 19951) and 
the Rayleigh limit again assumed. For silicate dust, corundum ab- 
sorption and scattering coefficients are also added to those of the 
silicates, with Jai = /si- This procedure is designed to approxi- 
mate the effect of heterogenous grains with a corundum core and 
a silicate mantle. Any pure corundum grains would have a much 
lower Td in the outer atmosphere than a grain that includes silicate 
with x < 1, and therefore would not be significant if hetrogenous 
grains are also present. 

Possibly, the optical constants of iKoike et alJ 1 19951) over- 
estimate the absorption in the near-infrared region of the spectrum, 
given that their optical constants k for amorphou s corun dum are 
high when compared to thin-film alumina iHarrisll 19551) or crys- 
talline corundum lHarman et aIlT 994). However, a lower value for 
k could not cause a > 20% change in the calculated stability ra- 
dius of corundum, as the gas is already becoming optically-thick in 
the mid-IR at the calculated stability radius due to molecular ab- 
sorption, meaning that the dust can not have a significantly lower 
temperature than the gas. Furthermore, a smaller stability radius 
for corundum would place it at radii smaller than those at which 
successful nucleation could occur (cf. Section l2~2l . 

2.4 Model Dust Types 

The dust types considered in this paper are given in TableQ Dust 
types A and B have a = 1.0, and represent olivine that will form 
pure forsterite initially in an outward-moving gas packet, before 
becoming enriched by Fe as the partial pressure of Mg decreases 
at / w 0.15. A higher value of a would give dust where Mg con- 
denses more fully before the addition of Fe to the olivine. The dif- 
ference between dust A and B is the number of grain nuclei, and 
hence the grain size: dust A will have dust grains a/TO times larger 
than dust B. Dust C is a dust where some the initial condensate 
contains an appreciable amount of Fe (x — 0.93), giving a higher 
absorption coefficient at wavelengths from 1 to 4/^m, where the 
bulk of the stellar flux is emitted. 

Dust type D, with corundum only, can be thought of as the 
dust that would exist if a < 0.4, as in this case silicates do not 
condense in our models. The effect of corundum on radiative trans- 
fer processes is at roughly the maximum plausible with the chosen 
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Table 1. Dust types considered in this paper. Columns are: Code used in this 
paper; Types of dust included (Olivine and Corundum); log of the number 
of grain nuclei per H nucleus; and a, the ratio of the exchange coefficient 
for Mg and Fe to the sticking coefficient for SiO (see text). 



Dust Code 


Dust Types 


log(AUc) 


a 


A 


Ol 


-13.3 


1.0 


B 


Ol 


-12.3 


1.0 


C 


Ol 


-12.3 


0.63 


D 


Cor 


-13 





value of 7V nuc = 10 , because much larger grains could not grow 
in an outward-moving gas packet even with sticking coefficients of 
1.0. 



3 MODELLING DETAILS AND RESULTS 

The atmospheric models used in this study are based on model se- 
ries of Hofmann et al. ( 1998), supplemented by models o et al. 
(2003b) and Ireland et al. 1 2004b c). These are non-grey dynamic 
model atmospheres for an M-type Mira variable with fundamental- 
mode period 332 days, mass 1.2 M® (M series; P series: 1.0 Mq) 
and solar metallicity. The luminosity L p and Rosseland radius R p 
of the non-pulsating parent star are 3470 Lp an d 260 R p, (M series; 
P series: 241 i? ). We refer to lHofmann et alj Il998h for essential 
details of model construction. In particular, these are self-excited 
models that were allowed to relax over a very large number of cy- 
cles. Pulsation is not strictly periodic, and the successive cycles 
that were studied in detail by Hofmann et al. (1998) and Ireland 
et al. (2004b, c) were so chosen that cycles of the M and P series 
with different characteristics are included. Phases assigned to these 
models (Ireland 2004b, c) are based on a reasonably chosen mean 
zero point, and the individual phase of a model of a specific cycle 
always refers to this zero point. The hydrodynamic models were 
calculated with dust-free grey opacities, and a non-grey tempera- 
ture stratification of the atmospheric layers was then obtained in 
a second step by solving the radiative-equilibrium equation (i.e. 
equation 1 with B v (T g ), gas absorption coefficient K vg (T g , P), 
gas temperature T g ) with dust-free non-grey opacities (see Hof- 
mann et al. 1998), retaining the gas-pressure stratification of the 
hydrodynamic model. 

During pulsation, the position of the continuum-forming lay- 
ers oscillates around Rp with an amplitude of ±30%. The param- 
eters for the models at the individual phases considered here are 
shown in Table|2| reproduced from Ire land et alj l2004bl) . For com- 
parisons with observations at specific visual phases, 0.12 will be 
added to (f> via in Table^l Irelan d et all2004bl) . 

Using these dust-free models as a starting point, the dust types 
as described in Section|2|were added to the models, and the tem- 
perature profile, spectra and intensity profiles re-calculated. The 
new gas-temperature T g and the dust-temperature Td stratifications 
were obtained by iteration, that is, radiative equilibrium was en- 
forced for both components (gas and dust) in each iteration step for 
obtaining for the next step new values of gas and dust temperatures, 
of gas and dust absorption coefficients , and of mean in tensities J„ 
(cf. comments on dust type D2 in Bedding et al. 2001). Note that 
scattering is approximated as being isotropic and coherent in these 
models and, hence, does not appear explicitly in Equation 1. Note 
also, that radiative acceleration does not become sufficiently large 
in the dusty models for generating a dust-driven wind, so that re- 



Table 2. Parameters of M series Mira models. The columns: visual phase 
(j> v i s \ luminosity L; 1.04 /im near-continuum radius Hi. 04; and the effec- 
tive temperature T\ 04 corresponding to iii.04- 



Model 


6 1 . 

"vis 


L 


^1.04 


Tl.04 






(Lq) 


(R P ) 


(K) 


1V1UJ 


fVi n AQ 


1470 


84 


2420 


Mf¥\n 
1V1UOI1 


fun £0 


2430 


78 


2860 


M08 


0+0 77 


4780 


81 


3320 


M09n 


0+0 89 


5060 


1 03 


2970 


M10 


1+0 02 


4910 


1 18 


2760 


ivi 1 1 n 


1 _i_n 1 1 








M12n 


1+0.21 


3470 


1.18 


2540 


M12 


1+0.27 


2990 


1.12 


2500 


M14n 


1+0.40 


1670 


0.91 


2400 


M15 


1+0.48 


1720 


0.83 


2530 


M16n 


1+0.60 


2460 


0.77 


2860 


M18 


1+0.75 


4840 


0.81 


3310 


M18n 


1+0.84 


4980 


1.00 


3010 


M19n 


1+0.90 


5070 


1.09 


2900 


M20 


2+0.05 


4550 


1.20 


2680 


M21n 


2+0.10 


4120 


1.21 


2610 


M22 


2+0.25 


2850 


1.10 


2490 


M23n 


2+0.30 


2350 


1.03 


2460 


M24n 


2+0.40 


1540 


0.87 


2410 


M25n 


2+0.50 


2250 


0.79 


2780 



1 0. 12 should be a dded to these phases to give a more accurate model phase. 
See llreland eTafl <2004bl) . 

taining the gas-pressure stratification of the original grey dynamical 
model is a reasonable approximation. 

The M series was chosen over the P series for this study be- 
cause approximations at 3-5 R p caused unrealistic density jumps 
with respect to time in the P series. These approximations included 
coarse gridding in the dynamical models, and an artificial density- 
gradient cutoff that had been enforced in the outer layers at some 
phases/cycles for computational reasons (non-physical density dis- 
continuity at the 5 R p model surface). The M series models consid- 
ered here have no artificial cutoff applied, and have a fine enough 
gridding of the dynamical models at 3-5 R p to produce no unreal- 
istic effects. 

The P and M model series are meant to describe typical Miras 
like o Cet or R Leo. Comparison with obs ervations o f both stars 
jHofmann et all200ltllreland et all2004jbHWoodruff et all2004 
Fedele et al. 2005J) show satisfactory agreement of observed and 
predicted feat ures, but various devia tions are also obvious from the 
discussion of Ireland et al. 1 2004c). We may check in this study 
whether some of these deviations are due to the omission of dust 
in the original models. Note that the mean-opacity treatment of 
strong water bands cools the outer layers too much and back-warms 
the con tinuum photosphere, as evident in the spectra of iTei et ail 
l2003bl) . If the forest of water lines were treated correctly in the 
models, one would expect the modelled water-bands to become less 
deep and the continuum (J-K) colour to become redder. This would 
give colour temperatures slightly lower than o Cet or R Leo. For 
this reason, we expect that realistic dust formation in these models 
would be at least as efficient as in o Cet or R Leo. 

FigureQshows the dust stability radii for all dust types in Ta- 
bleQ The dotted line shows the location of a single ma ss zone (the 
upper mass zone in Figure 2 of H ofmann et alj ll998l) ). giving an 
indication of the dynamics of the atmosphere. For dust types A, B 
and D, by far the largest rate of grain growth occurs just after initial 
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Figure 1. Model dust formation radii plotted against phase (0.12 has been 
added to the Table |2| phases - see text) for the four dust types considered 
here. The solid line shows the position where the dust-free near-continuum 
radius i?i.o4, and the dotted line shows the position of a single mass zone 
as a function of phase. 



grain formation, so it is in the initial stages of grain growth where it 
is most important to examine our assumption of dust condensation 
responding instantly to changes in the radiation field. For the mass 
zone shown in Figure □and zones immediately above it, the largest 
rate of grain growth occurs between phases 1.3 and 1.7, where the 
material crosses the dust stability radius. At these phases the mate- 
rial behind a strong shock (which reaches 2.75R P at phase 1.7) is 
cooling as it moves outwards. As the material crosses the dust sta- 
bility radius for dust types A and B at these phases, log(P) is be- 
tween — 1.9 and —1.7 with a local pressure scale height of 0.3R P . 
At these pressures, dust growth is sufficiently rapid according to 
the considerations of Section|2| However, this condition is not well- 
satisfied at all phases, so the details of the phase-dependence of the 
dust condensation radii should not be considered a reliable predic- 
tion of these simple models. 

The total optical depth of the dust at 1 .04 /im as a function 
of phase is shown in Figure [3] for all dust types. The large op- 
tical depths for dust type A are due to scattering by grains that 
reach a maximum radius of 96 nm in the outer layers, and reach 
80 nm radius rapidly after condensation while x is equal to 1 .0. The 
optical-depths for dust types A and B imply that the 1.04 /im near- 
continuum wavelength may not be a good window for interferomet- 
rically observing the continuum-forming photosphere as suggested 
by e.g. ljacob & Scholzl 120021) . As scattering opacity scales as A~ 4 
in the Rayleigh limit, the J and H band windows are not nearly as 
affected by the presence of these dust types as the 1 .04 /j,m window. 
It is clear that dust types C and D have would have only a minimal 
effect at this wavelength. 

The C-rich dust models o f H6f neretai] t2003). at nearly twice 
the luminosity of the M-series considered here, have a large radia- 
tive acceleration of the dust grains which significantly modifies the 
dynamic stratification of these models. At the 5 R p surface of the 
models considered here, the radiative acceleration ranges from 0.08 
to 0.29 times gravitational acceleration for dust A. This range of 
values is smaller for all other dust types or for radii less than 5 R p . 
This is easily understandable as dust only forms within 5 R p if the 
dust has a low opacity between 1 and 4 microns, where the bulk 
of the stellar flux is emitted. Therefore, we do not expect radiative 
acceleration of dust grains to play a highly significant role in the 
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Figure 2. Gas temperatures (lower lines) and dust temperatures (upper trun- 
cated lines) for the M22 models. The solid line is for the dust-free model, 
the dotted line dust A, the dashed line dust D and the dot-dashed line dust 
C. Dust type B is nearly indistinguishable from the dotted line and the gas 
temperature for dust C is nearly indistinguishable from the solid line. The 
dust-free near-continuum radius i?i.o4 for this model is 1.1 R p . 
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Figure 3. Total optical depths at 1 .04 urn as a function of phase for the four 
dust types. 



dynamic stratification of a typical Mira variable of luminosity less 
than about 5000Lq within 5 parent-star radii. However, beyond 
the surface of the models presented here radiative acceleration is 
expected to become significant as long as it is possible for the near- 
infrared opacity of the dust to increase with further Fe condensation 
at the low pressures (< 10 -4 dyn/cm 2 ) encountered at these radii. 

Figure [2] shows the influence of these dust types on the gas 
temperature for the M22 model. The gas in outer layers is warmed 
by about 100 K for dust types A and B. Back-warming isn't signifi- 
cant for any dust type in the layers where most spectral features are 
formed, which is at temperatures above about 1400 K. Other mod- 
els showed effects similar to those in this model. In this case, and 
in most other models, the dust temperatures are higher than the gas 
temperatures. This provides further justification for assuming that 
nucleation is complete when grain growth begins. The sharp change 
in slope of the dotted line at 2.5 R p is due to the the presence of 



6 M.J. Ireland and M. Scholz 



zero Fe (x = 1.0) for smaller radii, and an increasing fraction of 
Fe at larger radii. 



4 OBSERVABLE DUST EFFECTS 

The presence of all dust types considered here has a minimal effect 
on photometry between a wavelength of 1 and 9 /jm. Near mini- 
mum, the J band flux is reduced by 0.1 mag due to extinction by 
dust for dust types A and B, with no measurable effects on H and 
K band photometry. However, this small J-band effect is not a reli- 
able model prediction, as it requires growth and destruction rates in 
the outer layers near minimum that are too high for outer layer pres- 
sures P < 10 -3 dyn/cm 2 . Up to 0.3 mag extinction for dust type A 
occurs near-minimum at 900 nm due to back-scattered radiation re- 
intercepting the lower atmosphere. The extinction of V band light 
by the dust considered here is too difficult to model, both due to the 
difficulties in treating strong TiO absorption, and to the influence 
of the grain size distribution on extinction at short wavelengths. 

Interferometric observations at wavelengths short-wards of 
about 1 pm are sensitive to scattering from dust. Without allow- 
ing for significant effects of dust or molecules, lar ger than ex- 
pecte d apparent diameters have led some authors fe.g. lHaniff et alJ 
1995) to favour overtone pulsation for Mira variables. This pro- 
posa l is clearly in conflict with the now conve ntional understand- 
ing JWood et aljll999t IScholz & Woodll2000h that Miras pulsate 
in the fundamental m ode. Note that this d iscrepancy also existed 
in the near-infrared Ivan Belle et~ai]|l996t) . but correctly includ- 
ing the effects of molecular contamination of continuum diam- 
eters in ph ysical models has since provided a resolution in this 
case iWoodruff etaill2004 iFedele et alj feoOS). Several authors 
have clearly mentioned the possibility of dust scattering increas- 
ing the apparent size of Mira variables ( e.g. iDanchi et alJI 19941 : 
Hofm ann et"al]l200lt llreland et alj |2004d). but a physical model 
is required to discriminate correctly between the effects of TiO ab- 
sorption and dust. 

Table [3] shows measured diameters of o Cet and R Leo from 
several experiments compared with model predictions of the M22 
model with different dust types. The observation phases and tar- 
gets were: o Cet at phases 0.47 and 0.49 ( different cycles) fo r 
the William Herschel Telescope (WHT) data iTuthill et al1ll995h . 
o Cet a t phase 0.34 for the Anglo-Australian Telescope (AAT) data 
llreland et all2004 d1. R Leo at phase 0.2 for the Spe cial Astrophys- 
ical Observatory (SAO) data (Hofmann et al. 2001), and R Leo at 
phase 0.32 from Cambridge Optical Aperture Synthesis Telescope 
(COAST) lBurns etalJll998l) . The COAST data was presented as 
uniform-disk fits, and the correction factor of 1.52 in the text of 
iBurns etalJl 19981) was applied. 

For the model fits in Table [5] a least-squared Gaussian fit to 
visibility V (the normalised Fourier amplitude of the source bright- 
ness distribution) was calculated for baselines shorter than where 
V — 0.3. A Gaussian fit was both cited as a better fit than a uni- 
form disk by the authors of the observational papers, and is a bet- 
ter fit in general for the dusty models described here. The model 
star was placed at 105 pc, in-between the K-band maximum fit dis- 
tances for R Leo and o Cet for the M series llreland et alJ |2004b) 
and consistent with Hipparcos distances. Note that by having mod- 
els that match the observed temperature as measured by continuum 
J-K colours llreland et alj20 04b1 and using the K-band fit distance, 
the comparison between observed and predicted angular diameters 
is not heavily dependent on this distance. For example, if the mod- 
els are under-luminous by 20% with radii 10% too small, then the 



model stars will be placed at a distance 10% too close and the pre- 
dicted angular diameters will remain unchanged. 

For the 700, 750 and 920 nm filters, a 20 nm bandwidth was 
assumed, and for 1045 nm a monochromatic prediction is given. Al- 
though these bandwidths did not match the individual experimental 
bandwidths exactly, the experimental filters all included the chosen 
wavelengths and would be expec ted to have s imilar fit diameters 
according to the measurements of llreland et alJ l20 04d). 

It is clear that amongst these models, the observed large di- 
ameters and the increase in observed diameter at shorter wave- 
lengths can only be produced by dust A or B, with an intermediate 
dust being favorable. The scatter in the observations is indicative 
of both the slightly different phases of observations, and cycle-to- 
cycle variations in the atmospheric structure. This is also seen in 
the models, where the fitted FWHMs for the Ml 2 model would be 
smaller than those for the M22 model. 

The wavelengths of 700, 750, 920 and 1045 nm were chosen 
for this comparison because they are not sensitive to strong TiO ab- 
sorption in the upper atmosphere. In order to exclude TiO absorp- 
tion from causing the larger apparent diameters, we have carefully 
examined the treatment of TiO opacity in our models. The M series 
spectra have TiO absorption band depths that are much too deep 
when compared to obser vations. Although one significant reason 
for this in the figures of iTei et alJ J2003b) was out-dated chemi- 
cal da ta for Ti02 (since replaced with data from lSharp & Huebned 
1990), a remaining large problem in modelling the deep absorption 
bands correctly are non-LTE effects in TiO band formati on in a dy- 
namic atmosphere. From the TiO line list of lSchwenkd ll998l) . the 
typical Einstein A coefficient for the upper level in a strong TiO 
absorption band is 10 7 s -1 . This is much higher than the collision 
rate in the atmospheric regions where the TiO features are formed 
(P < 1 dyn/cm 2 and T < 2000 K). Several simple attempts have 
been made by us to characterize the magnitude of non-LTE effects. 
Indeed, spectra with band depths near that observed for typical Mi- 
ras can be produced by assuming an ad hoc TiO temperature profile. 
Although these tests affected the strong TiO features significantly, 
the regions of the spectrum such as those in Table [3] were largely 
un-affected. This is because these regions of weak and strongly 
temperature-dependent absorption are formed significantly deeper 
in the atmosphere than the strong absorption bands. Therefore, we 
are quite confident that the large apparent diameters can not be 
caused by TiO alone, and that therefore dust opacity is required 
to produce the observed diameters. 

The Opt i cal Int erferometric Polarimetry (OIP) observations of 
llreland et all 120051) are generally consistent with the predictions 
here. The radii of the dust shells in this paper are approximately 
2.5 R p , when the model star is placed at the distance corresponding 
to that which fits the K band maximum. Fitted optical-depths to 
these observations are between those predicted from the type A and 
type B dusts. However, the two stars observed in that paper (R Car 
and RR Sco) have smaller K-band photometric amplitudes than the 
M-series, and have much less dust emission than o Cet or R Leo. 
Therefore, one might expect the gas densities for these Miras to be 
smaller than that in the M series at the same radii. This would mean 
that grain formation and destruction may not be fast compared to 
the pulsation at the grain formation radii, and would need a more 
detailed, time-dependent study. 

In order to most accurately model the mid-infrared spectra, 
a new spectrum and intensity profile computation code was con- 
structed. This code used as inputs the gas pressure and molecule 
partial pressures, the gas and dust temperature and the velocity 
stratification from the model-construction code. It takes into ac- 
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Table 3. Model predictions for fitted Gaussian FWHMs to the stellar intensity profile in milli-arcsec for various wavelengths, for the M22 model which is at a 
phase of 0.37. Measurements for R Leo and o Cet for 4 separate instruments at visual phases between 0.2 and 0.49 are also shown (see text). 
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Figure 4. Model spectra for B dust between 8 and 13 microns for the Ml In 
model (lower solid line), M15 model (dotted line), Ml 8 model (dashed line) 
and M21n model (upper solid line). 



count continuum opaci ties, H2O lines fro m lpartridge & Schwenkd 
( 1997), TiO lines from Schwenke ( 1998) an d other dia tomic lines 
from the input to the ATLAS 12 models <KurucJfl994l) . The con- 
tributions to line profiles from different layers were redshifted or 
blueshifted as appropriate according to the difference in projected 
velocities. A micro-turbulence of 2.8km/s was assumed, in be- 
tween than us ed in the M-giant model s of IPlez et alj 1 19921) and 
that derived bv lHinkle & Barnes! <1979l) . 

The mid-infrared spectra are significantly influenced by dust 
types A and B, but the spectra for dust types C and D are similar 
to dust-free models. This is due to the low condensation fraction of 
dust C at radii where densities are high, and for dust D this is due 
to the low abundance of Al when compared to Si, Fe and Mg. As 
dust types A and B result in similar condensation fractions but dif- 
ferent grain sizes, their effects on the mid-infrared spectra are very 
similar. Hence, only selected mid-infrared spectra for representa- 
tive models with dust B are shown in Figure |4| smoothed to have 
a resolving power of 500. The model star is placed at a distance of 
102 pc, to match the fit distance for the o Cet K-band light curve 
faeland et all2004bl) . 

The most crucial point to note about these spectra is that they 
underestimate the measured flux from o Cet at all wavelengths 
by a factor of about 2 to 3, based on the range of photometry in 
Monni er et al. | <1998l) . This is consis ent with the measurements of 
Danc hi et al] 19941) and lWeiner et alj 1200 3l) . whose data require 
approximately half to three quarters of the total mid-infrared flux 
from o Cet to originate from a region further from the central star 
than our model surface (which corresponds to a 58 mas radius at 
a distance of 102 pc). This dust is almost certainly part of an out- 
flow from this star, and may have a shell-like structure as seen in 
outflows from other AGB stars (e.g. lHale et al Jfl997t) . Indeed, in 
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Figure 5. Visibility data used in Weine r et alj 120031) at a pulsation phase 
of 0.20 (data from November 28, 2000). One-parameter fits (see text) for 
models at similar phases (see text) are shown for the Ml In model (solid 
line) and the M21n model (dashed line). The dotted line illustrates a three- 
parameter fit (see text) based on the M21n model and an additional shell of 
previously ejected dust. 



the models a strong shock front occasionally (no more than about 
1 cycle in 4) drives a significant amount of mass through the 5 R p 
surface where it may form part of a stellar wind. Therefore, this 
additional flux is not inconsistent with the models here, it is simply 
not part of the models due to the arbitrary 5 R p model surface. 

A second important point to note about these model spectra is 
that they are not well approximated by a black-body spectrum and 
a silicate emission pe ak. According to the classification scheme of 
ISloan & Pricd fr998). the spectra for the M series with dust B could 
be classified from SE2 to SE8, with a dust emission contrast of 0. 14 
to 0.3. This suggests that the observed mid-infrared excess for some 
Mira variables may not be strictly associated with an outflow (e.g. 
Ilreland et all2005l) . 

The effect of these dust types on mid-infrared interferometry 
is difficult to predict. This is due to the large fraction of the light 
that originates fro m dust further than the 5 R p surface to these mod- 
els iWeiner et al]l2003t) . Nevertheless, one can say that a general 
feature of the models is that they do not resemble a uniform-disk 
intensity profile. Figure|5| shows 1 1.15 fim V 2 data from o Cet ob- 
servations at a phase of 0.2 obtained from J. Weiner (personal com- 
munication ). This is calibrated data which was used for the model 
fits in IWeiner et alj J5)03), unlike the un-calibrated data plotted 
in Figure 1 of that paper. Intensity profiles were computed using 
dust B, and the narrow observing bandwidth, assuming a redshift 
of 83km/s. Over-plotted in Figure [5] are the model predictions of 
the Ml In model, at a pulsation phase of 0.23, and the predictions 
of the M21n model at a pulsation phase of 0.22. The model star 
was placed at the fit distance for the o Cet K-band light curve 
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Jlrelandetall ExHb) and the fraction of emission from an over- 
resolved component of extended dust emission was a free parame- 
ter. This fraction has the effect of multiplying the squared visibility 
by a constant smaller than 1. 

Certainly in the case of the Ml In model the fit is reasonable, 
and the 45 % additional flux required from extended dust emission 
is roughly consistent with observed mid-infrared excesses for o Cet. 
In the case of the M21 model, additional model parameters would 
be required to fit the data reasonably. The example plotted as the 
dotted line is one where two additional parameters were free: the ra- 
dius of an additional shell of dust and the fractional emission from 
this shell. The parameters for this model were a dust shell radius of 
5.83 R p containing 18 % of the total emission. This demonstrates 
that for mid-infrared interferometry to provide accurate constraints 
on the innermost dust and water emission, a model or observations 
of dust at large radii must be combined with models of the inner re- 
gions as described here. Examining data from several other epochs 
did not give significant errors when fitting dust B models to the 
data, but dust-fr ee models were not in general extended enough to 
fit the data from Wei ner et alJ<2003l) . 

A final observable feature of the dust formation models as de- 
scribed here is the phase- and cycle-dependence of the dust for- 
mation. A relatively quick increase in dust opacity is predicted 
near minimum during certain cycles wherever a shock front passes 
through the radius of dust formation. This occurs around phase 1.4 
in these models, where the 1.04 /im dust optical depth rapidly in- 
creases from 0.26 to 0.72 for dust A. This increase in opacity would 
be best seen interferometrically by a jump in size at either wave- 
lengths short-wards of 1 /im or in the mid-infrared. Combined ob- 
servations at these two wavelength regimes could constrain grain 
sizes (e.g. the difference between dust types A and B). 

The cycle-dependence of the outer layers in the M and P se- 
ries relates to their chaotic nature. There is no reason to expect these 
layers to show high degrees of spherical symmetry as assumed in 
the models. A crude approximation is that one could think of op- 
posite sides of a star to be in different model 'cycl es'. Sim ilar be- 
haviour is evident on large scales in the models of Woifkd l2005t) . 
For this reason, dust formation near-minimum may also be observ- 
able as the creation of asymmetric features in a shell at 2-3 R p . 
One effect of this dust creation would be rapid changes in the angle 
of polarization observed. This effect is a plausible explanation for 
the rapid changes in polarization angl e observed a round Julian day 
2441000 (near phase 0.5) for R Leo bv lSerkowski & Shawlli200ll) . 
These observations show polarizations of several percent in U, B 
and V bands where the angle of polarization changes by 25 degrees 
in one direction over 50 days, then 50 degrees in the opposite di- 
rection over the following 50 days. 



5 CONCLUSION 

Condensation of corundum and silicate dust in a model of a 1 .2 Mq 
self-excited Mira variable has been examined. Four dust types were 
modelled, representing approximated physical dust formation pa- 
rameters. It was found that a dust type between A and B from Ta- 
ble Q fits existing observations well . In particular, the effect of 
scattering by this type of dust explains the systematically large ob- 
served apparent diameters at wavelengths shorter than 1 /im. This 

1 Predictions of existing and additional models are available upon request 
to anyone interested for specific observational programs. 



dust type has approximately 10~ grain nuclei per H atom avail- 
able at the onset of silicate condensation, and requires a, the ratio 
of the exchange coefficient for Mg and Fe ion to the sticking coef- 
ficient of SiO, to be at least the order of unity. A model where only 
corundum forms was found to be a poor fit to observations of o Cet 
and R Leo. The use of composition-dependent opacities enabled 
the survival of the initial Mg-rich silicate condensate at radii of 2-3 
times the model parent star radius. 
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